Interplay of soft and hard processes and hadron p T spectra in pA and AA collisions 
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Motivated by a schematic model of multiple parton scattering within the Glauber formalism, the 
transverse momentum spectra in pA and AA collisions are analyzed in terms of a nuclear modification 
factor with respect to pp collisions. The existing data at the CERN Super Proton Synchrotron (SPS) 
energies are shown to be consistent with the picture of Glauber multiple scattering in which the 
interplay between soft and hard processes and the effect of absorptive processes lead to nontrivial 
nuclear modification of the particle spectra. Relative to the additive model of incoherent hard 
scattering, the spectra are enhanced at large pr (hard) by multiple scattering while suppressed at 
low pr (soft) by absorptive correction with the transition occurring at around a scale po ~ 1 — 2 
GeV/c that separates soft and hard processes. Around the same scale, the pr spectra in pp collisions 
also change from an exponential form at low px to a power-law behavior at high px- At very large 
Pt > Po, the nuclear enhancement is shown to decrease like l/Px- Implications of these nuclear 
effects on the study of jet quenching, parton thermalization and collective radial flow in high-energy 
A A collisions are discussed. 
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I. INTRODUCTION 



Hadron yields, spectra and correlations have been the 
focus of many experiments in relativistic heavy-ion col- 
lisions. They provide a snapshot of the state of mat- 
ter when hadrons stopped interacting with each other, 
a stage often referred to as freeze-out in heavy-ion col- 
lisions. One can then infer the condition in the early 
stage prior to the freeze-out. Such a procedure relies 
crucially on our understanding of the dynamical evolu- 
tion of the system. On the other hand, knowledge of the 
initial condition at the beginning of thermalization will 
also help one to unravel the history of evolution. This is 
especially the case when a complete thermalization can- 
not be achieved in certain regions of phase space (large 
px, for example). Therefore, it is important to study the 
hadron spectra in p + p and p+A collisions that will help 
us to understand the initial condition of the dense matter 
in high-energy heavy-ion collisions. 

Hadron production at large transverse momentum in 
p + p collisions is relatively well understood within the 
perturbative QCD (pQCD) parton model. Because of 
the large transverse momentum scale, the hard parton- 
parton scattering processes can be calculated in pQCD, 
while the non-perturbative effects can be factorized into 
universal parton distributions and parton fragmentation 
functions. These parton distributions and fragmenta- 
tion functions can be independently measured in deep- 
inelastic electron-nucleon scattering, e + e~ annihilation 
and other processes. The calculated inclusive hadron 
spectra either in leading order jjj or next-to-leading order 
H agree very well with experimental data. 

In pA collisions, one should take into account multiple 
scatterings inside the nucleus. Because of the interfer- 
ence effect, these multiple scatterings will give rise to 
different A-scaling behavior of the spectra at different 
values of px- Generalizing to AA collisions, one should 



expect similar behavior for the spectra of initially pro- 
duced partons. These produced partons would undergo 
further interactions, possibly leading to thermalization 
(or partial thermalization) and expansion, which would 
be reflected in the final hadron spectra. Therefore, one 
should understand first the nuclear modification of the 
hadron spectra from p + p to p + A and to A + A due 
to initial multiple scattering. Only then can one reliably 
extract information about the freeze-out conditions, e.g., 
temperature and flow velocity and parton energy 

loss H from the final hadron spectra. 

In this paper, we will analyze the nuclear modification 
of the hadron pt spectra in p + A collisions motivated 
by a schematic model of multiple parton scattering in 
which coherence and absorptive corrections are included 
via the Glauber multiple scattering formalism. Multiple 
scatterings generally enhance the large px spectra rel- 
ative to the additive model of hard scattering. It can 
be shown || that absorptive corrections and the power- 
law behavior of perturbative parton cross section are the 
main reasons why the nuclear enhancement at high pt 
decreases like l/p T - We will also show that the same 
absorptive processes suppress the spectra relative to the 
additive model at low px, where soft processes dominate 
and the px spectra deviate from a power-law behavior. 
By analyzing the experimental data in terms of the pro- 
posed nuclear modification factor, we can phenomenolog- 
ically determine the px scale that separates soft and hard 
processes by the transition from nuclear suppression to 
nuclear enhancement of the px spectra. 



II. MULTIPLE PARTON SCATTERING IN P + A 
COLLISIONS 



Multiple hard parton scatterings in QCD are consid- 
ered as high-twist processes and their contributions to 
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the cross section of hadronic collisions are generally sup- 
pressed by 1/Q 2 , where Q is the momentum scale of the 
processes. The coefficients of the these contributions are 
generally related to the multiple parton correlation func- 
tions inside a hadron In collisions involving a nu- 
cleus, such high-twist contributions are enhanced by A 1 / 3 
due to the large nuclear size. For certain quantities like 
the transverse momentum imbalance of dijets in photon- 
nucleus collisions || , parton energy loss in deep- inelastic 
lepton-nucleus scattering ||, broadening of the jet trans- 
verse momentum in deep- inelastic lepton-nucleus scatter- 
ing and Drell-Yan pairs in p + A collisions |ll[] , these 
high-twist terms are the leading contributions. However, 
such QCD treatment of multiple parton collisions is so far 
limited to the parton level and has not yet been extended 
to hadron spectra. In this paper, we will use instead a 
schematic Glauber model 12 1 of multiple parton scatter- 
ing to motivate our analysis. 

Let us denote h 3 iN as the differential cross section for 
a parton-nucleon scattering, i + N — ► j + X, 
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where f b /N{ x b) is the parton distribution function inside 
a target nucleon, s, t and u are Mandelstam variables 
and d&i b ^j c /dt are the differential parton-parton cross 
sections. In terms of h 3 iN , we define the effective parton- 
nucleon total cross section as 



(2) 



and the differential nuclcon-nucleon cross section for par- 
ton production as 
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The effective parton-nucleon and nucleon-nucleon cross 
sections are only finite after one introduces some effec- 
tive infrared cut-off in the parton-parton scattering pro- 
cesses. We will discuss this cut-off later. Following the 
approach by [fl2|]6f| , in which the Glauber approximation 
is used for multiple parton scattering, one can find the 
parton spectra in pA collisions up to double scattering 
approximation as 
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where Ja(&) is the nuclear thickness function with nor- 
malization J d 2 bTji{b) — A. We take the nucleus to be 
a hard sphere of radius ta — r^A 1 / 3 with ro = 1.14 fm. 



One can then obtain the ratio of the differential cross 
sections oiN + A^j + X and N + N -> j + X, 
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Compared to the additive model of parton scattering 
which gives Ra = 1, the second term in the above equa- 
tion gives us the nuclear modification of the parton spec- 
trum due to multiple parton scattering inside a nucleus. 
This term contains contributions from both the double 
parton scattering and the absorptive correction to the 
single scattering processes. Notice that the absorptive 
contribution at this order is negative. As we will demon- 
strate, it is the cancellation by the absorptive correction 
that leads to many interesting and nontrivial features of 
the nuclear modification of the particle spectra. 



III. NUCLEAR MODIFICATION: 
STUDY 



A SCHEMATIC 



As a demonstration of the consequences of multiple 
parton scattering, absorptive correction and the resul- 
tant interplay between low and high pt behavior of the 
particle spectra in pA collisions, we illustrate the nuclear 
modification of parton spectra in a schematic model. As- 
suming parton-hadron duality, the conclusions can be ap- 
plied qualitatively to hadron spectra. 

For a schematic study, let us assume that all partons 
are identical and the differential parton-nucleon cross sec- 
tions have a simple regularized power-law form in px, 



h J lN = h{p T ) = 
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(P 2 T 



Pi) 



(\y\ < AF/2), 



with the total parton-nucleon cross section, 
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Such a form is motivated by the fact that both the pQCD 
calculation and experimentally measured differential jet 
production cross sections show such a power-law behav- 
ior. The parameter pq is introduced as an infrared cut- 
off. It can be regarded phenomenologically as a scale that 
separates soft and hard processes, because for p T < p 
the parton spectra deviate significantly from a power-law 
behavior. 
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Using this schematic parton-nucleon cross section, the 
nuclear modification factor can be simplified as 
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This modification factor as a function of Pt/po f° r differ- 
ent values of n is shown in Fig. [l], where we plot i?^ — 1 in 
units of the modification strength < dA 1 / 3 a / '(167TT 2 ,). One 
can also evaluate the nuclear modification factor analyt- 
ically at two different limits of pr, 
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FIG. 1. The nuclear modification factor as a function of 
Pt for different values of n in a schematic model of multiple 
parton scattering with a simple form of parton-nucleon cross 
section l/[pT +Po) n - 

As one can see from both Fig. ^ and the above asymp- 
totic behavior, the parton spectra are enhanced at large 
Pt due to contribution from double scattering. However, 
due to the cancellation by the absorptive correction to the 
single scattering, the enhancement decreases like Pq/p t 
at large pt- This asymptotic behavior will never occur 
if one does not include the absorptive corrections which 
are embedded in the Glauber formalism employed in this 
model. At low Pt, contributions from double scattering 
become smaller than the absorptive correction, so that 



the effective parton spectra are reduced (Ra < 1) rela- 
tive to the additive model. For very large values of the 
pp cross section and the size of the nucleus A, we notice 
that the modification factor Ra could become negative in 
Eq. ([)]) at pt — 0, signaling the breakdown of the double 
scattering approximation. In this case, one would have 
to include triple and other multiple scatterings. One can 
show along the same line as the double scattering, the 
contribution from the sum of triple scattering and the 
corresponding absorptive correction is positive at pt = 
and is proportional to (A 1 / 3 ^) 2 . Assuming one can expo- 
nentiate the sum of all multiple scatterings and absorp- 
tive corrections, the modification factor at pt = would 
take the form, 
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In this regime, the absorptive parts dominate and Ra oc 
1/A 1 / 3 . This is consistent with the wounded nucleon 
model for soft particle production. The coefficient c < 1 
is a result of partial cancellation by contributions from 
multiple scattering. 

It is clear that the shape of the nuclear modification of 
the inclusive spectra results from the interplay between 
multiple scattering and absorptive corrections. This is 
especially true at low pt- Even though it gives the ap- 
parent broadening of the effective spectra in pA relative 
to pp collisions, the underlying physics is very different 
from the random-walk model |13j for the pt broadening 
of soft hadrons. One can expect the same behavior in 
A + A collisions if there is no additional nuclear effect, 
such as final state rescattering or jet quenching due to ra- 
diative parton energy loss. Analyzing the experimental 
data in this framework with respect to the above base- 
line behavior would allow one to identify and study new 
effects caused by the presence of dense matter. 

Another important feature of the modification factor 
is its sensitivity to the form of particle spectra in pp col- 
lisions. At low pt, its behavior is not very sensitive to 
the form of differential parton-nucleon cross section. If 
a Gaussian form of the spectrum is usded, the modifica- 
tion factor remains roughly the same at pt ~ 0. In con- 
trast, the high-pT behavior of the enhancement, ~ l/p T , 
is strictly the consequence of the power-law form of the 
parton-nucleon cross section. In fact, with a Gaussian 
form, the enhancement will increase with p\ exponen- 
tially. Therefore, the observed large pt behavior of the 
Cronin enhancement JTifl is consistent with the power- 
law form of the jet production cross section and with 
the picture of multiple parton scattering in p + A colli- 
sions. The nuclear modification factor at large pt also 
depends on the power n of the pp spectra. Since the 
power n decreases with energy, one should expect Ra to 
decrease with energy. This trend has been observed in 
experiments 15 in the energy range y/s =20-40 GeV. 

Overall the nuclear modification factor Ra has nontriv- 
ial and interesting pt dependence. It is smaller than 1 at 
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low pt and larger than 1 at intermediate pr- The tran- 
sition (Ra = 1) occurs at around pt ~ Po/V™ ~ VvPr)- 
It approaches 1 again at very large pr Po- The modi- 
fication strength at both low and high pt is proportional 
to the size of the nucleus. As we will see in the next 
section, this is qualitatively consistent with the experi- 
mental data of inclusive hadron spectra in pA collisions. 
Since po is a momentum scale that set the onset of power- 
law-like spectra in pp collisions, one can consider it a 
scale separating soft and hard processes. According to 
the schematic model, one should be able to determine 
independently this scale from the nuclear modification 
factor of the hadron spectra. 



IV. ANALYSIS OF EXPERIMENTAL DATA 

To facilitate the analysis of pA and AA spectra in terms 
of the nuclear modification factor Ra, one needs to know 
the pp spectra at the same energy. Since we will have 
to take ratios of hadron spectra from different experi- 
ments with different pt bins, this can only be achieved 
by parameterizing the baseline spectra. Shown in Fig. ^ 
are the measured negative hadron spectra |16 17j] and the 
parameterization in pp collisions at E\ a \> = 200 GeV. We 
use a two-component parameterization, 



f( PT ) = C e- m ^ To + C 



(1 - 2 PT /^sr 



(P 2 T 
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The fit parameters are listed in Table | for different spec- 
tra at different energies. The parameters in the expo- 
nential form are mainly determined by the spectra below 
Pt < 1 GeV/c while those in the power-law form are 
mainly determined by the large pt region. Shown as 
a dashed line is the contribution from the exponential 



component. It is clear that at pr > 2 GeV/c, the spec- 
trum is already dominated by the power-law behavior. 
The power-law form we use for the CERN-SPS energy 
range contains a factor (1 — Ipxj \fs) a that is caused by 
the rapid decrease of quark distributions at large x ~ 1. 
This factor will become negligible at higher energies when 
2pr/V^ < 1. 
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FIG. 2. Negative hadron spectra in pp collisions at 
Blab = 200 GeV from Refs. and the fit (solid line) ac- 

cording to Eq. ( |ll| ) with fit parameters given in Table |l[ The 
dashed line is the underlying exponential component. The 
hg is defined as = (ir + + vr~)/2 + K~ + p. The upper 
inserted box shows the ratio of parameterizations (solid line) 
at Slab = 158 and 200 GeV and the corresponding pQCD 
parton model calculation (dot-dashed line). The lower box 
shows the ratio of h~ /■k and h$ /ir . 




TABLE I. Fit parameters of the hadron spectra in pp collisions 





Co (mbGeV" 2 ) 


To (MeV) 


C (rnbGeV" 2 ) 


a 


po (GeV/c) 


n 


h-(pp,E uh = 158 GeV) 


169.6 


154 


24565 


9.3 


2.43 


6.29 


h-{pp,E lab = 200 GeV) 


174.7 


154 


13653 


9.1 


2.27 


6.22 


h-(pp 1 E Ub = WQ GeV) 


174.7 


154 


17653 


9.8 


2.27 


6.29 


7v^(pp,E lRb = 200 GeV) 


150.3 


154 


13653 


9.1 


2.37 


6.22 


^(pp, ^8 = 200 GeV) 


440 


154 


653 





1.75 


4.98 



In the original UA1 parameterization llq], a power-law 
form A 1 (pt +Po) n is used to fit the hadron spectra over 
a large range of energies. While this single power-law 
is sufficient to fit the spectra for pt > 0.2 GeV/c, an 
exponential term is necessary to fit the spectra at low 
transverse momentum px < 0.2 GeV/c. 

In order to compare to the spectra in pA and AA colli- 
sions for near isospin symmetrical nuclei, the data shown 
in Fig. H at large pr are for " negative" hadrons with the 



7r~ contribution being replaced by an isospin averaged 
value, Zip = (tt + +tt~)/2 + K~ +p. Since particle produc- 
tion at large px is dominated by the leading hadrons from 
valence quark scattering and there are more up-quarks 
than down-quarks in a pp system, one should see more 
7r + than 7r~ at large pt- Shown in the inserted box at 
the lower-left corner are the ratios of h~ /n^ and /n^ 
[7r ± = (tt + +7t~)/2]. Note that h~ /7r ± decreases with pt 
as expected while /n^ remains relatively constant. At 
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low pt, h~ and are approximately the same within 
a few percent accuracy. This isospin dependence can be 
described well by the pQCD parton model |^3| . With Zip 
and h~ , one will be able to estimate the isospin effect in 
pA and AA collisions. In the inserted box at the upper- 
right corner we also plot the ratio of the hadron spectra 
in pp collisions at E\ a b = 158 and 200 GeV. Since there 
are no experimental data for pp collisions at Eiab = 158 
GeV, we parameterize the spectrum at i?i a b = 158 GeV 
to fit the ratio (dot-dashed line) obtained from a pQCD 
parton model calculation p3| . 
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and the fit (solid line) accord- 



FIG. 3 
yfs = 200 GeV from Ref. |_ 



ing to Eq. ([n]) with fit parameters given in Table |l[ The 
dashed line is the underlying exponential component and the 
dot-dashed line is the power-law component. 

At the BNL Relativistic Heavy-ion Collider (RHIC) 
energy, we expect the power-law contribution to become 
more important. Shown in Fig. |^ is the charged hadron 
spectrum in pp collisions at *Js = 200 GeV |1J] and the 
parameterization. The power-law is indeed more promi- 
nent than at the CERN-SPS energy. It becomes domi- 
nant already at around pt — 1-5 GeV/c. There is a gen- 
eral trend that both po and n of the power-law component 
decrease with colliding energy. However, the underlying 
exponential term remains the same. 

With these parameterizations of hadron spectra in pp 
collisions, we can analyze the hadron spectra in pA and 
AA collisions in terms of a nuclear modification factor 
which is defined in general as 



Rab{pi 



d(JAB/dyd 2 p T 



(N hina , Ty }d(TNN/dyd 2 pj 



(12) 



for AB collisions, where (Aginary) = / d 2 bTAB(b) is the 
number of binary collisions averaged over the impact- 
parameter range of the corresponding centrality. Here 



TabQ>) is the nuclear overlap function for AB collisions 
at impact parameter b. For minimum-biased events, 
(Abinary) = AB . For the purpose of the study in this 
paper, we will select event centrality according to the 
geometrical cross sections. 
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FIG. 4. The nuclear modification factor Rab{pt) for 
hadron spectra in pA and AB collisions at _Ei a b = 200 GeV. 
The data are from Refs. 100,01. 



Shown in Figs. | and | are the nuclear modification 
factors for hadron spectra in pA and AB collisions with 
different centralities at CERN-SPS energies. One can see 
that both pA and AB data have features very similar to 
those shown by the schematic model. The nuclear mod- 
ification factors are smaller than 1 at low px- According 
to Eq. (|l0|), the modification factor at pt = 0, Rab(0), 
should decrease with increasing nuclear size. This is in- 
deed seen in the data. As pt increases, Rab also in- 
creases and becomes larger than 1 at about pr = 1 — 2 
GeV/c. According to the physical picture depicted by the 
schematic model, this momentum scale represents the on- 
set of hard processes underlying the hadron spectra. The 
increase of Rab is then due to the onset of the hard com- 
ponent of incoherent parton scattering underlying the 
hadron production mechanism as pt increases. The data 
then indicate that at pt = 1 — 2 GeV/c, the hard compo- 
nent becomes completely dominant. Indeed, such a value 
is consistent with the momentum scale around which the 
hadron spectra in pp collisions begin to be dominated 
by a power-law behavior of hard parton scatterings. In 
principle, the nuclear modification factor should peak at 
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intermediate values of px and decrease again to approach 
to 1 at very large px- However, at pt = 4 GeV/c, one is 
already close to the kinetic limit of y/s/2 at the CERN- 
SPS energy. Near the kinetic limit the spectra become 
sensitive to other nuclear effects such as the Fermi mo- 
tion that will increase the parton distribution function at 
x r-j 1, leading to increase of Rab again. In pA collisions 
at the RHIC energy, one will certainly see the decrease 
of Ra again at large px ^> 2 GeV/c. 



10 



b 

A 



V 1 

\ 

CD 

+ 

b 



E lob =158 AGeV Pb + Pb 

A tt° 10% central (WA98) 
• h" 5% central (NA49) 
A 7t 8% central Pb+Au (CERES) 




10 1 

Pt (GeV/c) 

FIG. 5. The nuclear modification factor Rab(pt) for 
hadrons in AB collisions at Eiat = 158 GeV. The data are 
from Refs. pp|-p2[. The line is the pQCD parton model cal- 
culaiton. 



The increase of the nuclear modification factor with 
Pt due to the onset of hard parton scatterings appears 
very similar to the broadening of hadron spectra due to 
collective radial flow in a hydrodynamic picture. If there 
are strong interactions among partons in the early stage 
and hadrons in the late stage of heavy-ion collisions, 
the system will be driven to local thermal equilibrium. 
The information of initial multiple parton scattering con- 
tained in the initial parton spectra could be partially or 
completely erased by the thermalization. Collective flow 
then will be developed and hadron spectra will become 
broader than in pp and pA collisions due to the boost by 
a collective flow velocity [ p4| . Such an effect essentially 
will also increase Rab as a function of px- The question 
is how one can distinguish these two apparently different 
dynamics that produce the same final hadron spectra. 
Mass dependence of the hadron spectra has been pro- 
posed as a unique measure of the collective flow effect in 
heavy- ion collisions J2^]. Since a collective flow provides 
a common velocity boost for all particles, heavy parti- 
cles will then acquire more transverse momentum in the 
nonrelativistic region (px ^ rrih). One then should see a 
linear mass dependence of the slope parameter from an 
exponential fit of the measured hadron spectra pq |. In 
the current parton model p3|] , one cannot exclude the 
mass dependence of the nuclear modification factor. It is 
therefore important to perform the current analysis for 



different hadron species in pp and pA collisions to find 
out whether nuclear modification factors have a mass de- 
pendence. Only then can one find out whether there is 
complete or partial thermalization and to what extent 
the effect of initial multiple parton scatterings has sur- 
vived the final state interaction and contributed to the 
apparent collective radial flow. 



V. HIGH P T SPECTRA AND PARTON ENERGY 
LOSS 

Since the interaction for an energetic parton inside a 
medium is dominated by small angle forward scatterings, 
thermalization will be less complete for high px partons. 
Hydrodynamic description will then become less relevant. 
Large px hadron spectra will be determined by how an 
energetic parton propagates inside the medium. Since 
the large px hadron spectra are calculable in the pQCD 
parton model, they are good probes of the parton dy- 
namics in dense matter. The nuclear modification factor 
at large px is a convenient and efficient way to deter- 
mine the effects of the final state interaction of energetic 
partons inside a medium. 
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FIG. 6. The predicted nuclear modification factor Rab(jpt) 
for 7ro spectra in Pb + Pb collisions at _Ei a b = 158 GeV with 
different centrality cuts. 

To demonstrate whether the experimental data at the 
CERN-SPS are consistent with the picture of multiple 
parton scattering, we also plot in Fig. [| the pQCD- 
inspircd parton model calculation of the nuclear modi- 
fication factor In such a model, one extends the 
collinear factorized parton model to include intrinsic 
transverse momentum and its broadening due to mul- 
tiple scattering in nuclear matter. The value of the in- 
trinsic transverse momentum and its nuclear broadening 
are adjusted once and the model can reproduce most of 
the experimental data in pp and pA collisions p3| . In AA 
collisions, the effect of parton energy loss is modeled by 
the modification of parton fragmentation functions |2q] 
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in which the distribution of leading hadrons from parton 
fragmentation is suppressed due to parton energy loss. 
We will not describe the model here and refer readers to 
Ref. HjJ for detail. We note that the behavior of the nu- 
clear modification factor at px > 2.0 GeV/c in Pb + Pb 
collisions is well described by the model without any ad- 
ditional final state medium effect such as the parton en- 
ergy loss or jet quenching ||. Since the pQCD parton 
model cannot deal with soft processes, the calculations 
are only reliable for px > 2 GeV/c. The px broadening 
due to multiple parton scattering is essential to account 
for the pt dependence of the nuclear modification factor. 
Different variations of the parton model studies ]27j give 
the same conclusion. Such a good agreement is not just 
a coincidence. The calculated absolute differential cross 
section for different colliding systems (S + S, S + Au 
and Pb + Pb) and for different centralities also agree well 
with the experimental data p3fl . Shown in Fig. ^ are the 
calculated nuclear modification factors for 7r° spectra in 
Pb + Pb collisions with different centrality cuts, assum- 
ing no parton energy loss. Though the px dependences 
of the modification factors are similar, the magnitudes 
increase with centrality. The impact-parameter ranges 
for different centralities are determined by the fractions 
of geometrical nuclear cross sections. 

Recent theoretical studies of parton propagation in a 
dense medium predict a substantial parton energy loss 
due to induced gluon radiation as the parton interacts 
with the dense medium ]28|-|32]] . If a dense medium is pro- 
duced in heavy-ion collisions, large px jets will lose energy 
as they propagate through the medium. This will lead to 
suppression of large px hadrons However, as shown 
by Fig. the experimental data of central Pb + Pb col- 
lisions at the CERN-SPS are consistent with the picture 
of initial multiple parton scattering without any parton 
energy loss in the final state. Careful analysis of the data 
against the pQCD parton model can actually exclude ef- 
fects of any significant amount of parton energy loss in 
central Pb + Pb collisions at the CERN-SPS ||. To be 
consistent with the current theoretical estimates of par- 
ton energy loss in a dense medium, these data imply that 
the initial energy density would be small and the lifetime 
of the dense system is very short in Pb + Pb collisions 
at the CERN-SPS. On the other hand, preliminary data 
from RHIC experiments |m| have shown significant sup- 
pression of high px hadron spectra. Combined with the 
null results on hadron suppression in heavy-ion collisions 
at the CERN-SPS, RHIC data would imply the onset of 
parton energy loss due to the increasing initial energy 
density and the lifetime of the dense system created in 
heavy-ion collisions at RHIC energies. Study of the nu- 
clear modification factors for hadron spectra at large px 
will help to determine the average parton energy loss in 
the medium produced in heavy- ion collisions |3| . 

The nuclear modification factor will also help us to 
understand the energy dependence of the parton energy 
loss. Earlier theoretical studies [p8|-|30"[ gave a constant 
or a weakly energy-dependent parton energy loss dE/dx. 



However, more recent studies in an opacity expansion ap- 
proach |n],[52| predict a stronger energy-dependence for 
dE/dx when the parton energy is small (E < 10//, \x 
being the Debye screening mass in the medium). Pa- 
rameterization of the numerical results of Ref. gives 
dE/dx cx ^ 15 /(10 + E) 1 - 25 for fx = 0.5 GeV. The coef- 
ficient depends linearly on the initial parton density and 
the system size. This form has a strong energy depen- 
dence for E ~ 10 GeV and a weak dependence at large 
E. 



Au+Au(b = 0) Vs=130GeV 




Pt (GeV/c) 

FIG. 7. The parton model calculation of nuclear modifica- 
tion factor Rab{pt) for hadrons in central Au + Au collisions 
at y^s = 130 GeV with different forms of parton energy loss. 
A q is the value of mean free path for a propagating quarks 
used in the calculation. 

Shown in Fig. fj] are the nuclear modification factors 
calculated in a pQCD parton model for different 
forms of parton energy loss. The nuclear modification 
factor is suppressed to become smaller than 1 for non- 
vanishing parton energy loss. For a constant dE/dx, 
the modification factor increases with px after the ini- 
tial drop. This is because that the constant energy loss 
becomes relatively less important for higher initial par- 
ton energy. For asymptotically large px, the finite and 
constant energy loss would become negligible and the 
modification factor Rab would approach 1. However, 
for an energy loss that has a strong energy dependence, 
the shape of the modification factor is very different. If 
this is the case, the energy loss for a parton in the small 
px region is still very small and only becomes sizable at 
large px- Consequently, the effect of energy loss is also 
very small. The modification factor will then increase, 
following the trend of the initial parton spectra due to 
the transition from soft to hard scatterings. However, 
when larger energy loss sets in at higher px, the hadron 
spectra are strongly suppressed leading to much smaller 
values of Rab- We should note that the dE/dx used in 
these calculations is only representative of its value aver- 
aged over the entire evolution of the dense system. Since 
the energy loss is directly proportional to the parton den- 
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sity of the medium which decreases very fast due to rapid 
expansion, a small averaged dE/dx still corresponds to 
large initial energy loss and large initial parton density. 

Similar to the interplay between soft and hard pro- 
cesses in the initial parton production, there should also 
be a smooth transition between effects of hydrodynamic 
evolution at low px and parton energy loss at high px- 
The hydrodynamic picture is relevant only if there is lo- 
cal thermalization, which is more likely for low px par- 
tons. It breaks down for very large px partons because 
the small angle scatterings in QCD are not sufficient to 
bring these partons in equilibrium with the rest of the 
system. These high px partons will suffer energy loss 
as they propagate through the medium leading to sup- 
pression of large px hadrons. However, at intermediate 
px, one cannot neglect the effect of partial thermaliza- 
tion or the detailed balance. For example, one should 
consider both gluon absorption and radiation by a prop- 
agating parton. So far all theoretical studies have only 
considered gluon radiation. If gluon absorption is also 
included, the effective parton energy loss will be reduced 
|{35f . Since the effect of gluon absorption will be small 
for high energy partons (px 3> T, T being the tempera- 
ture of the thermal medium) , the absorption will further 
increase the energy dependence of the effective parton 
energy loss. The px dependence of the nuclear modifica- 
tion factor Rab(pt) a t intermediate px will then provide 
useful information about the thermalization of the dense 
medium. 



VI. CONCLUSIONS 

In this paper we have proposed to analyze the hadron 
transverse momentum spectra in terms of the nuclear 
modification factor Rab{pt) which is defined in such a 
way that a naive additive model of incoherent hard par- 
ton scattering would give Rab — 1- We demonstrated 
in a schematic model of Glauber multiple parton scatter- 
ing that the modification factor Rab{pt) has a nontriv- 
ial px dependence due to the absorptive processes and 
the interplay between soft and hard parton scattering, 
excluding final state scatterings. Because of the absorp- 
tive processes, the hadron production at small px ~ 
is coherent and the hadron spectra in AB collisions are 
proportional to the number of participant nucleons, lead- 
ing to Rab < 1. At large px the hard parton scatterings 
become incoherent. Multiple parton scatterings then en- 
hance the hadron spectra so that Rab > 1- The momen- 
tum scale po at which the transition occurs [Rab{po) = 1] 
can be identified as the scale that separates soft and hard 
processes underlying both pp and AB collisions. Analy- 
ses of the existing experimental data on pp, pA and AB 
collisions indicate that po w 1 — 2 GeV/c. 

We pointed out that such analyses of future experi- 
mental data are important to study the effect of final 
state interactions. At low px, collective radial flow from 



hydrodynamic expansion gives similar px dependence of 
Rab ■ Disentangling the effects of initial multiple scatter- 
ing and the radial flow would require a careful study of 
the modification factor in pA collisions, especially its de- 
pendence on the hadron mass. At large px, parton energy 
loss will lead to suppression of the hadron spectra. Ex- 
perimental measurement of Rab will provide important 
information on the initial parton density that is produced 
in heavy- ion collisions. At the intermediate px, we have 
shown that the px dependence of Rab is sensitive to the 
energy dependence of dE/dx. This in turn is related to 
gluon absorption by the propagating partons reflecting 
the detailed balance in an equilibrating system. 

Although experimental data at the CERN-SPS H do 
not indicate any sign of parton energy loss, recent pre- 
liminary experimental data from RHIC |34| show signifi- 
cant suppression of large px hadrons. Detailed studies of 
the nuclear modification factor Rab{pt) over the whole 
range of px in both pA and A A collisions, will help us 
not only to determine the parton energy loss or the ini- 
tial parton density but also to find out the degree of 
thermalization and radial collective flow. Recent stud- 
ies Q also pointed out that parton energy loss can also 
lead to azimuthal anisotropy in high px hadron spectra 
which is significantly different from hydrodynamic behav- 
ior. A combination of studies on Rab and the azimuthal 
anisotropy at high px will provide a window to the initial 
condition and dynamics of early evolution of the dense 
matter that was never possible before. 
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